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EXISTENCE OF INDEPENDENT RANDOM MATCHING 

By Darrell Duffie and Yeneng Sun 

Stanford University and National University of Singapore 

This paper shows the existence of independent random match- 
ing of a large (continuum) population in both static and dynamic 
systems, which has been popular in the economics and genetics liter- 
atures. We construct a joint agent-probability space, and randomized 
mutation, partial matching and match-induced type-changing func- 
tions that satisfy appropriate independence conditions. The proofs 
are achieved via nonstandard analysis. The proof for the dynamic 
setting relies on a new Fubini-type theorem for an infinite product 
of Loeb transition probabilities, based on which a continuum of in- 
dependent Markov chains is derived from random mutation, random 
partial matching and random type changing. 

1. Introduction. Economists and geneticists, among others, have implic- 
itly or explicitly assumed the exact law of large numbers for independent 
random matching in a continuum population, by which we mean a nonatomic 
measure space of agents. This result is relied upon in large literatures within 
general equilibrium theory, game theory, monetary theory, labor economics, 
illiquid financial markets and biology, as discussed in [15], which provides 
extensive references. Such a law of large numbers allows a dramatic sim- 
plification of the dynamics of the cross-sectional distribution of properties 
among a large population of randomly interacting agents. Mathematical 
foundations, however, have been lacking, as has been noted by Green and 
Zhou [19]. 

Given the fundamental measurability problems associated with modeling 
a continuum of independent random variables, 1 there has, up to now, been 
no theoretical treatment of the exact law of large numbers for independent 
random matching among a continuum population. In [40], various versions 
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of the exact law of large numbers and their converses are proved by apply- 
ing simple measure-theoretic methods to an extension of the usual prod- 
uct probability space that has the Fubini property. 2 The measure-theoretic 
framework of [40] is adopted in our companion paper [15] to obtain an ex- 
act law of large numbers for random pairwise matching by formulating a 
suitable independence condition in types. 3 

In particular, assuming independent random matching in a continuum 
population, and some related independence assumptions regarding random 
mutation and match-induced type changing, we prove in [15] that there 
is an almost sure deterministic cross-sectional distribution of types in a 
large population for both static and dynamic systems, a property that had 
been widely used without a proper foundation. In addition, we show in [15] 
that the time evolution of the cross-sectional distribution of types can be 
completely determined from the agent-level Markov chain for type, with 
known transition matrices. 

The main aim of this paper is to provide the first theoretical treatment of 
the existence of independent random matching in a continuum population. 
In particular, we construct a joint agent-probability space, and randomized 
mutation, partial matching and match-induced type-changing functions that 
satisfy the independence conditions in [15]. Though the existence results of 
this paper are stated using common measure-theoretic terms, their proofs 
make extensive use of nonstandard analysis. One can pick up some back- 
ground knowledge on nonstandard analysis from the first three chapters of 
the book [34]. 

Since our dynamic system of independent random matching generates a 
dependence structure across time (in particular, a continuum of indepen- 
dent Markov chains of agents' type), we need to construct an infinite prod- 
uct of Loeb transition probabilities. Our proof of the existence result for 
the dynamic setting is based on such an infinite product and an associated 
Fubini- type theorem that is derived from Keisler [26]. Specifically, in order 



2 These results were originally stated on Loeb measure spaces in [38]. However, as noted 
in [40] , they can be proved for an extension of the usual product space that has the Fubini 
property; see also Chapter 7 in [34] (and in particular. Sections 7.5 and 7.6), written by 
Sun. 

3 The independence condition that we propose in [15] is natural, but may not be obvious. 
For example, random matching in a finite population may not allow independence among 
agents since the matching of agent i to agent j implies of course that j is also matched 
to i, implying some correlation among agents. The effect of this correlation is reduced to 
zero in a continuum population. In other words, one may obtain standard independence by 
rounding infinitesimals but not *-independence for a hyperfinite population (see Section 4 
below). A new concept, "Markov conditional independence in types," is proposed in [15] 
for dynamic matching, under which the transition law at each randomization step depends 
on only the previous one or two steps of randomization. 
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to prove the Fubini-type property of an infinite product of Loeb transition 
probabilities, we first generalize Keisler's Fubini theorem for the product of 
two Loeb probability measures, from [26], to the setting of a Loeb transition 
probability (i.e., a class of Loeb measures) in Theorem 5.1. In construct- 
ing the internal transition probabilities for the step of random mutation in 
Section 6, we use a hyperfinite product space and its coordinate functions, 
which are closely related to those of Keisler [26], and in particular to the 
law of large numbers for a hyperfinite sequence of *-independent random 
variables, as noted in [26], page 56. 

Historically, reliance on the exact law of large numbers for independent 
random matching dates back at least to 1908, when G. H. Hardy [20] and W. 
Weinberg (see [6]) independently proposed that random mating in a large 
population leads to constant and easily calculated fractions of each allele in 
the population. Hardy wrote: "suppose that the numbers are fairly large, 
so that the mating may be regarded as random," and then used, in effect, 
an exact law of large numbers for random matching to deduce his results. 
Consider, for illustration, a continuum population of gametes consisting of 
two alleles, A and B, in initial proportions p and q = 1 — p. Then, following 
the Hardy- Weinberg approach, the new population would have a fraction p 2 
whose parents are both of type A, a fraction q 2 whose parents are both of 
type B, and a fraction 2pq whose parents are of mixed type (heterozygotes). 
These genotypic proportions asserted by Hardy and Weinberg are already, 
implicitly, based on an exact law of large numbers for independent random 
matching in a large population. We provide a suitable existence framework. 

Going from a static to a dynamic environment, we also provide an exis- 
tence result that allows the computation of a steady-state constant deter- 
ministic population distribution of types. For illustration, suppose in the 
Hardy- Weinberg setting above that with both parents of allele A, the off- 
spring are of allele A, and that with both parents of allele B, the offspring are 
of allele B. Suppose that the offspring of parents of different alleles are, say, 
equally likely to be of allele A or allele B. The Hardy- Weinberg equilibrium 
for this special case is a population with steady-state constant proportions 
p = 60% of allele A and q = 40% of allele B. With the law of large numbers 
for independent random matching, this is verified by checking that, if gen- 
eration k has this cross-sectional distribution, then the fraction of allele A 
in generation k + 1 is almost surely 0.6 2 + 0.5 x (2 x 0.6 x 0.4) = 0.6. Our 
existence results for a dynamic model of random matching provide a math- 
ematical foundation for this Hardy-Weinberg law governing steady-state 
allelic and genotypic frequencies. 

In the field of economics, Hellwig [21] is the first, to our knowledge, to have 
relied on the effect of the exact law of large numbers for random pairwise 
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matching in a market, in a 1976 study of a monetary exchange economy. 4 
Since the 1970s, a large economics literature has routinely relied on an exact 
law of large numbers for independent random matching in a continuum 
population. This implicit use of this result occurs in general equilibrium 
theory (e.g., [17, 18, 35, 41]), game theory (e.g., [4, 5, 8, 16, 22]), monetary 
theory (e.g., [11, 19, 21, 31]), labor economics (e.g., [10, 23, 36, 37]) and 
financial market theory (e.g., [14, 32]). Mathematical foundations, however, 
have been lacking, as has been noted by Green and Zhou [19]. In almost all 
of this literature, dynamics are crucial. For example, in the monetary and 
finance literature cited above, each agent in the economy solves a dynamic 
programming problem that is based in part on the conjectured dynamics 
of the cross-sectional distribution of agent types. An equilibrium has the 
property that the combined effect of individually optimal dynamic behavior 
is consistent with the conjectured population dynamics. In order to simplify 
the analysis, much of the literature relies on equilibria with a stationary 
distribution of agent types. 

The remainder of the paper is organized as follows. Section 2 consid- 
ers existence of independent random matching, both full and partial, in a 
static setting, after a brief introduction of the measure-theoretic framework 
in Section 2.1. Theorem 2.4 of Section 2.2 shows the existence of random 
full matching with independent types, meaning roughly that, for essentially 
every pair of agents, the type of the agent to be randomly matched 
with agent i is independent of the type of the agent to be randomly matched 
with agent j. Theorem 2.6 of Section 2.3 then considers existence for the 
case of random search and matching, that is, for random partial matchings 
that are independent in types. Proofs of Theorems 2.4 and 2.6, which use 
nonstandard analysis extensively in the computations, are given in Section 
4. 

Section 3 considers a dynamical system for agent types, allowing for ran- 
dom mutation, partial matching and match-induced random type changes. 
We borrow from our companion paper [15] the inductive definition of such 
a dynamical system given in Section 3.1, and the condition of Markov con- 
ditional independence found in Section 3.2. The latter condition captures 
the idea that at every time period, there are three stages: (1) an indepen- 
dent random mutation, (2) an independent random partial matching, and 



4 Diamond [9] had earlier treated random matching of a large population with, in effect, 
finitely many employers, but not pairwise matching within a large population. The match- 
ing of a large population with a finite population can be treated directly by the exact law 
of large numbers for a continuum of independent random variables. For example, let N(i) 
be the event that worker i is matched with an employer of a given type, and suppose this 
event is pairwise independent and of the same probability p, in a continuum population 
of such workers. Then, under the conditions of [40] , the fraction of the population that is 
matched to this type of employer is p, almost surely. 
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(3) for those agents matched, an independent random type change induced 
by matching. In economics applications, for example, match-induced type 
changes arise from productivity shocks, changes in asset holdings induced by 
trade between the matched agents, changes in credit positions, or changes in 
money holdings. Theorem 3.1 of Section 3.3 shows the existence of a dynam- 
ical system B with random mutation, partial matching and type changing 
that is Markov conditionally independent in types with any given param- 
eters. Theorem 3 of [15] then implies that the type processes of individual 
agents in such a dynamical system B form a continuum of independent 
Markov chains, and that the time evolution of the cross-sectional distribu- 
tion of types is deterministic and completely determined from a Markov 
chain with explicitly calculated transition matrices. 5 

We prove Theorems 2.4 and 2.6 in Section 4. Turning to Section 5, we 
first prove in Section 5.1 a generalized Fubini theorem for a Loeb transi- 
tion probability. Then, a generalized Ionescu-Tulcea theorem for an infinite 
sequence of Loeb transition probabilities is shown in Section 5.2. Finally, 
a Fubini extension based on Loeb product transition probability systems 
is constructed in Section 5.3. Based on the results in Section 5, we prove 
Theorem 3.1 in Section 6. 

Finally, we emphasize again that we must work with extensions of the 
usual product measure spaces (of agents and states of the world), since a 
process formed by a continuum of independent random variables is never 
measurable with respect to the completion of the usual product c-algebra, 
except in the trivial case that almost all of the random variables in the 
process are constants. 6 

2. Existence of independent random matchings in the static case. In this 
section, we first give some background definitions in Section 2.1. Then, we 
consider the existence of random matchings that are independent in types, 
for full and partial matchings, in Sections 2.2 and 2.3, respectively. Proofs 
of the two existence theorems, Theorems 2.4 and 2.6, are given in Sections 
4.1 and 4.2. The full and partial matching models in Theorems 2.4 and 2.6 
satisfy the respective conditions in Theorems 1 and 2 in [15], which means 
that the respective conclusions of Theorems 1 and 2 in [15], characterizing 
the implications of the law of large numbers for process describing the cross- 
sectional distribution of types, hold for these matching models. 



5 The models in [27, 28] and [29] assume that there is a set of individual agents and 
a single marketmaker. At each time, a randomly chosen individual agent trades with the 
marketmaker. References [1] and [2] formalize a link between matching and informational 
constraints, which do not consider random matching under the independence assumption 
as in the model here. 

6 See, for example, Proposition 1.1 in [39]. 
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2.1. Some background definitions. Let probability spaces (1,1, X) and 
(QjTjP) be our index and sample spaces, respectively. 7 In our applications, 
(1,1, A) is an atomless probability space that indexes the agents. 8 Let (/ x 
Q,T®T, X®P) be the usual product probability space. For a function / on 
/ x ft (not necessarily J (g) ^-measurable) , and for (i, u) & I x£l, represents 
the function f(i,-) on ft, and the function f(-,u>) on /. 

In order to work with independent type processes arising from random 
matching, we need to work with an extension of the usual measure-theoretic 
product that retains the Fubini property. A formal definition, as in [40], is 
as follows. 

Definition 2.1. A probability space (/ x f2, W,Q) extending the usual 
product space (/ x J 7 , X®P) is said to be a Fubini extension of (/ x O, 

T®T , \ <£>P) if for any real- valued Q-integrable function g on (/ x £2, W), the 
functions gi = g(i, ■) and g u = /(•, uj) are integrable respectively on (fi, J 7 , P) 
for A-almost all i £ i and on (1,1, A) for P-almost all u> £ Q; and if, moreover, 
j n g-i dP and J 7 g u dX are integrable respectively on (I,T, A) and on (f2,.F, P), 
with J IX Qg dQ = Jj(J n gidP) dX = J n (fj g^dX) dP. To reflect the fact that 
the probability space (/xSl, W, Q) has (I, I, A) and (Q, T , P) as its marginal 
spaces, as required by the Fubini property, it will be denoted by (/ x Vt,ZM T , 
A KIP). 

An X M ^"-measurable function / will also be called a process, while ft is 
called a random variable of the process and is called a sample function 
of the process. 

We now introduce the following crucial independence condition. We state 
the definition using a complete separable metric space X for the sake of 
generality; in particular, a finite space or a Euclidean space is a complete 
separable metric space. 

Definition 2.2. An Z Kl JT-measurable process / from / x to a com- 
plete separable metric space X is said to be essentially pairwise independent 
if for A-almost all i £ I, the random variables fi and fj are independent for 
A-almost all j € /. 



7 A11 measures in this paper are countably additive set functions denned on cr-algebras. 

8 A probability space (1,1, A) is atomless if there does not exist AgT such that X(A) > 
0, and for any X-measurable subset C of A, A(C) = or A(C) = X(A). For those interested 
in the case of a literal continuum of agents, it is noted in the beginning of Section 4 that 
one can indeed take / to be the unit interval with some atomless probability measure. 
Corollary 4.3 in [40] shows, however, that, in general, it makes no sense to impose the 
Lebesgue measure structure when an independent process is considered. 
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2.2. Existence of independent random full matchings. We follow the no- 
tation in Section 2.1. Below is a formal definition of random full matching. 

Definition 2.3 (Pull matching). 

1. Let S = {1,2,..., K} be a finite set of types, a: I — > S an T-measurable 
type function of agents and p its distribution on S. For 1 < k < K, let 
Ik = {i £ I : a(i) = k} and p^ = A(2fc) for each 1 < k < K. 

2. A full matching <p is a bijection from I to I such that for each i £ I, 
(/>(i) 7^ i and cfr((j)(i)) = i. 

3. A random full matching ir is a mapping from / x O to / such that (i) tt w 
is a full matching for each u £ f2; (ii) if we let 5 be the type process a(7r), 
then 5 is measurable from (/ x Q,I MT, A IE P) to S 1 ; (hi) for A-almost 
all i £ I, gi has distribution p. 

4. A random full matching 7r is said to be independent in types if the type 
process g is essentially pairwise independent. 

Condition (1) of this definition says that a fraction p^ of the population 
is of type k. Condition (2) says that there is no self-matching, and that if 
i is matched to j = (j>(i), then j is matched to i. Condition (3) (Hi) means 
that for almost every agent i, the probability that i is matched to a type-A; 
agent is pk, the fraction of type-A: agents in the population. Condition (4) 
says that for almost all agents i and j € /, the event that agent i is matched 
to a type-A; agent is independent of the event that agent j is matched to a 
type-? agent for any k,l £ S. 

The following theorem shows the existence of an independent random full 
matching model that satisfies a few strong conditions that are specified in 
footnote 4 of [35], and is universal in the sense that it does not depend on 
particular type functions. 9 Note that condition (1) (ii) below implies that for 
any i, j £ I, P{fti = j) = since A({j'}) = 0, which means that the probability 
that agent i is matched with a given agent j is zero. 

Theorem 2.4. There exists an atomless probability space (1,2, X) of 
agents, a sample probability space (£l,J-,P), a Fubini extension (I x $l,ZM 
J-,\E\P) of the usual product probability space, and a random full matching 
vr from (I x n, 1M T, A M P) to I such that 



9 When (1,1, A) is taken to be the unit interval with the Borel algebra and Lebesgue 
measure, property (l)(ih) of Theorem 2.4 can be restated as "for P-almost all lo G fi, 
A(Ai HttJ 1 ^)) = \(Ax)\(Ai) holds for any A\ , A2 £ X" by using the fact that the count- 
able collection of rational intervals in [0, 1] generates the Borel algebra. Footnote 4 of [35] 
shows the nonexistence of a random full matching it that satisfies (i)-(iii) of part (1). 
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1. (i) for each lo £ $7, X(7r UJ 1 (A)) = X(A) for any A el, (ii) for each i € 
/, P{n~ l (A)) = \(A) for any A el, (iii) for any Ai,A 2 € 1, \{A X n 
tt~ 1 (A 2 )) = \{A{]\{Ai) /"5^(is /or JF '-almost all lo G $7; 

2. 7r is independent in types with respect to any given type function a from 
I to any finite type space S. 

2.3. The existence of independent random partial matchings. We shall 
now consider the case of random partial matchings. The following is a formal 
definition. 



Definition 2.5. Let a: J— > S be an X- measurable type function with 
type distribution p = (pi, . . . ,px) on S. Let it be a mapping from I x Q to 
I L) {J}, where J denotes "no match." 

1. We say that tt is a random partial matching with no-match probabilities 
qi,...,q K in [0, 1] if: 

(i) For each w£fl, the restriction of n w to / — vrJ 1 ({J}) is a full 
matching on / — 7rJ 1 ({J}). 10 

(ii) After extending the type function a to ID { J} such that a( J) = J, 
and letting g = a(-7r), we have g measurable from (/ x £1,1 M J 7 , A M P) to 
SU{J}. 

(iii) For A-almost all i £ Ik, P(gi = J) = Qk and 11 

Er=lPr(l-9r) 

2. A random partial matching 7r is said to be independent in types if the 
process g (taking values in SU { J}) is essentially pairwise independent. 12 

The following theorem generalizes Theorem 2.4 to the case of random par- 
tial matchings. Because the given parameters for no-matching probabilities 
may be type-dependent, it is not possible to produce a universal matching 
model for random partial matchings as in the case of full matchings. 



10 This means that an agent i with 7r w (i) = J is not matched, while any agent in / — 
7rJ ({J}) is matched. This produces a partial matching on /. 

11 If an agent of type k is matched, its probability of being matched to a type-? 
agent should be proportional to the type distribution of matched agents. The fraction 
of the population of matched agents among the total population is ^2^= 1 Pr(l — q r )- 
Thus, the relative fraction of type-? matched agents to that of all matched agents is 
— Qi))/ ^2^-iPr(^ ~ Qr)- This implies that the probability that a type-fc agent is 

matched to a type-Z agent is (l — q k )(pi(l — qi))/Y^=iP^0-~1^)- When E^liP'^ 1 _ Qr) = 
0, we have — qk) = for all 1 < k < K, in which case almost no agents are matched, 
and we can interpret the ratio ((1 — qk)pi0-~ gi))/y^ J L 1 Pr(l — Qr) as zero. 

12 This means that for almost all agents i,j £ /, whether agent i is unmatched or 
matched to a type-fc agent is independent of a similar event for agent j. 
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Theorem 2.6. There is an atomless probability space (1,1, A) of agents 
such that for any given 1 -measurable type function (3 from I to S, and 
for any q S [0, l] s , (1) there exists a sample space (£l,T,P) and a Fubini 
extension (I x Q,Z M T, A M P) of the usual product probability space; (2) 
there exists an independent-in-types random partial matching ir from (I x 
0,XKUF, A M P) to I with q = (qi, . . . ,qx) as the no-match probabilities. 

3. The existence of a dynamical system with random mutation, partial 
matching and type changing that is Markov conditionally independent in 
types. A discrete-time dynamical system B with random mutation, partial 
matching and type changing that is Markov conditionally independent in 
types is introduced in Section 4 of our companion paper [15]. The purpose 
of this section is to show the existence of such a dynamical system ID with 
any given parameters. In Sections 3.1 and 3.2, we reproduce respectively 
the inductive definition of a dynamical system with random mutation, par- 
tial matching and type changing and the condition of Markov conditional 
independence, which originated with [15]. The general existence of the dy- 
namical system B is presented in Theorem 3.1 in Section 3.3 and its proof 
in Section 6. 

3.1. Definition of a dynamical system with random mutation, partial match- 
ing and type changing. Let S = {1,2, ... ,K} be a finite set of types. A 
discrete-time dynamical system B with random mutation, partial matching 
and type changing in each period can be defined intuitively as follows. The 
initial distribution of types is p°. That is, p°(k) (denoted by p£) is the initial 
fraction of agents of type k. In each time period, each agent of type k first 
goes through a stage of random mutation, becoming an agent of type I with 
probability 6jy. In models such as [14], for example, this mutation generates 
new motives for trade. Then, each agent of type k is either not matched, 
with probability qk, or is matched to a type-Z agent with a probability pro- 
portional to the fraction of type-Z agents in the population immediately after 
the random mutation step. When an agent is not matched, she keeps her 
type. Otherwise, when a pair of agents with respective types k and I are 
matched, each of the two agents changes types; the type-k agent becomes 
type r with probability Vki( r ), where uti is a probability distribution on S, 
and similarly for the type-/ agent. 

We shall now define formally a dynamical system B with random muta- 
tion, partial matching and type changing. As in Section 2, let (1,1, A) be 
an atomless probability space representing the space of agents, (Q,,T,P) a 
sample probability space, and (I x £1,1 M T, A M P) a Fubini extension of 
the usual product probability space. 

Let a : 1 — *■ S = {1, . . . , K} be an initial T-measurable type function with 
distribution p° on S. For each time period n > 1, we first have a random 
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mutation that is modeled by a process h n from (I x £1,1 M J 7 , XM P) to 
S, then a random partial matching described by a function 7r n from (/ x 
£1,2 M T , A M P) to I L) {J} (where J represents no matching), followed by 
type changing for the matched agents that is modeled by a process a n from 
(I x Q,ISf,A§P) to S. 

For the random mutation step at time n, given a K x K probability 
transition matrix 13 b, we require that, for each agent i £ I, 

(1) P(h» = l\ c %- 1 = k)=b kl , 

the specified probability with which an agent i of type k at the end of time 
period n — 1 mutates to type I. 

For the random partial matching step at time n, we let p™- 1 / 2 be the ex- 
pected cross-sectional type distribution immediately after random mutation. 
That is, 

(2) plr l/2 =p n - l/2 {k)= J^X({ieI:K(i)=k})dP(u). 

The random partial matching function n n at time n is defined by: 

1. For any lu £ £1, 7r"(-) is a full matching on / — (vr2)~ 1 ({ J}), as defined in 
Section 2.3. 

2. Extending /i n so that h n (J,uj) = J for any w£fi, we define g n :I x £1^ 
SU{J} by 

ff R (t,u;) = /i n (7r n (i J w) ) a;) ) 

and assume that g n is X (El ^"-measurable. 

3. Let g £ [0, 1] . For each agent i £ J, 

W = W = *0 = <fc. 

(3) 

p( g r=^r^) = (1 7 fc , )(1 " 9 5vT - 

Equation (3) means that for any agent whose type before the matching 
is k, the probability of being unmatched is qk, and the probability of being 
matched to a type-Z agent is proportional to the expected cross-sectional type 
distribution for matched agents. When g n is essentially pairwise independent 
(as under the Markov conditional independence condition in Section 3.2 
below), the exact law of large numbers in [38] and [40] implies that the 
realized cross-sectional type distribution A(/i") _1 after random mutation at 
time n is indeed the expected distribution p n-1 / 2 , P-almost surely. 14 



14 



Here, bki is in [0, 1], with bu = 1 for each k. 

As noted in Footnote 11, if the denominator in (3) is zero, then almost no agents will 



be matched and we can simply interpret the ratio as zero. 
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Finally, for the step of random type changing for matched agents at time 
n, a given v. S x S ^ A specifies the probability distribution v^i = v{k, I) of 
the new type of a type-Zc agent who has met a type-£ agent. When agent i is 
not matched at time n, she keeps her type hf with probability 1. We thus 
require that the type function a n after matching satisfies, for each agent 
i € I, 

P(a? = r\h? = k,g? = J)=6l 

(4) 

P(a? = r\h? = k,g? = l) = v kl (r), 

where 5£ is 1 if r = k, and zero otherwise. 

Thus, we have inductively defined a dynamical system D with random mu- 
tation, partial matching and match-induced type changing with parameters 
(p°,b,q,u). 

3.2. Markov conditional independence in types. In this subsection, we 
consider a suitable independence condition on the dynamical system B. For 
n > 1, to formalize the intuitive idea that given their type function a n_1 , 
the agents randomly mutate to other types independently at time n, and 
that their types in earlier periods have no effect on this mutation, we say 
that the random mutation is Markov conditionally independent in types if, 
for A-almost all i G I and A-almost all j € /, 

P(fc? = k, h] = l\al a?" 1 ; «?, • • • , a]- 1 ) 

(5) 

= P(h? = k\a™- l )P{h™ = lla]- 1 ) 

holds for all types k, I £ S 1 . 15 

Intuitively, the random partial matching at time n should depend only on 
agents' types immediately after the random mutation. One may also want 
the random partial matching to be independent across agents, given events 
that occurred in the first n — 1 time periods and the random mutation at 
time n. We say that the random partial matching 7r n is Markov conditionally 
independent in types if, for A-almost all i € / and A-almost all j El, 

P(gf = c,g? = d\al a?" 1 , hf; a°, . . . , a]~\h]) 

(6) 

= P(g? = c\h?)P(g? = d\h]) 

holds for all types c, d G S U { J}. 

The agents' types at the end of time period n should depend on the 
agents' types immediately after the random mutation stage at time n, as 



15 We could include the functions h m and g m for 1 < m < n — 1 as well. However, it is 
not necessary to do so since we only care about the dependence structure across time for 
the type functions at the end of each time period. 
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well as the results of random partial matching at time n, but not otherwise 
on events that occurred in previous periods. This motivates the following 
definition. The random type changing after partial matching at time n is 
said to be Markov conditionally independent in types if for A-almost all i £ I 
and A-almost all j G I, and for each n > 1, 

) • • • > i^i i9i i ®j > ■ • ■ ) &j > hj , 5j ) 

holds for all types k,l E S. 

The dynamical system D is said to be Markov conditionally independent 
in types if, in each time period n, each random step (random mutation, 
partial matching and type changing) is so. 

3.3. The existence theorem. The following theorem shows the existence 
of a dynamical system with random mutation, partial matching and type 
changing that is Markov conditionally independent in types. Its proof will 
be given in Section 6 after the development of a generalized Fubini theorem 
for a Loeb transition probability and a Loeb product transition probability 
system in Section 5. 

Theorem 3.1. Fixing any parameters p° for the initial cross-sectional 
type distribution, b for mutation probabilities, q £ [0, l] s for no-match prob- 
abilities and v for match-induced type-change probabilities, there exists a 
Fubini extension of the usual product probability space on which is defined a 
dynamical system B with random mutation, partial matching and type chang- 
ing that is Markov conditionally independent in types with these parameters 
(p°,b,q,u). 

Note that the dynamic matching model D described in Theorem 3.1 above 
satisfies the conditions in Theorem 3 of [15]. Thus, the conclusions of Theo- 
rem 3 in [15] also hold for the matching model D, including the statements 
that the type processes of individual agents in such a dynamical system D 
form a continuum of independent Markov chains, and that the time evo- 
lution of the cross-sectional type process is completely determined from a 
Markov chain with known transition matrices. 

4. Proofs of Theorems 2.4 and 2.6. In this section, we give the proofs 
for Theorems 2.4 and 2.6 in Sections 4.1 and 4.2, respectively. 

The space of agents used from this section onward will be based on a 
hyperfinite Loeb counting probability space (1,1, X) that is the Loeb space 
(see [33] and [34]) of internal probability space (I,Io,\q), where I is a hy- 
perfinite set, Iq its internal power set, and Aq(^4) = |A|/|7| for any A^Iq 
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(i.e., Ao is the internal counting probability measure on I). Using the usual 
ultrapower construction as in [34], the hyperfinite set / itself can be viewed 
as an equivalence class of a sequence of finite sets whose sizes go to infinity, 
and the external cardinality of / is the cardinality of the continuum. Thus, 
one can also take the unit interval [0, 1] as the space of agents, endowed with 
a (j-algebra and an atomless probability measure via a bijection between / 
and the unit interval. 

All of the internal probability spaces to be discussed from this section on- 
ward are hyperfinite internal probability spaces. A general hyperfinite inter- 
nal probability space is an ordered triple ($7, JF , Po), where Q = {u\, u-z, . . . ,lo^} 
for some unlimited hyperfinite natural number 7, is the internal power 
set on $7, and Pq{B) = Yli<j<y,u . e # PoiWj}) f° r an Y B € Fq. When the 
weights Po({ w j})> 1 < j < 7) are all infinitesimals, (£1,J-q,Pq) is said to be 
atomless, and its Loeb space (fi,^ 7 , P), as a standard probability space, is 
atomless in the usual sense of Footnote 8. Note that nonstandard analysis 
is used extensively from this section onward. The reader is referred to the 
first three chapters of [34] for more details. 

4.1. Proof of Theorem 2.4. Fix an even hyperfinite natural number in 
the set *Noo of unlimited hyperfinite natural numbers. Let I = {1, 2, . . . , N}, 
let Iq be the collection of all the internal subsets of /, and let Ao be the 
internal counting probability measure on Iq. Let (1,1, X) be the Loeb space 
of the internal probability space (I,2o,Xq). Note that (1,1, X) is obviously 
atomless. 

We can draw agents from / in pairs without replacement, and then match 
them in these pairs. The procedure can be the following. Take one fixed 
agent; this agent can be matched with N — 1 different agents. After the first 
pair is matched, there are N — 2 agents. We can do the same thing to match 
a second pair with N — 3 possibilities. Continue this procedure to produce a 
total number of 1 x 3 x • • • x (N - 3) x (N - 1), denoted by (N - 1)!!, different 
matchings. Let 17 be the space of all such matchings, Tq the collection of 
all internal subsets of $7 and Po the internal counting probability measure 
on Tq. Let (Q,,!F,P) be the Loeb space of the internal probability space 

(n,r ,p Q ). 

Let (I x Q,1q <g> Tq, Ao <8> -Po) be the internal product probability space of 
(I,Tq, Ao) and (Q,,Tq,Pq). Then Zq ® is actually the collection of all the 
internal subsets of / x Q and Ao <8> -Po is the internal counting probability 
measure on Iq® J-q. Let (/ x T , A KIP) be the Loeb space of the internal 

product (I x Q,Io 8) T§, Ao ® Po), which is indeed a Fubini extension of the 
usual product probability space. 16 



16 For any given two Loeb spaces (1,1, X) and (Sl,F,P), it is shown in [30] that the 
Loeb product space (I x f2,XEl T, A M P) is uniquely defined by the marginal Loeb spaces. 



14 



D. DUFFIE AND Y. SUN 



Now, for a given matching u G SI and a given agent i, let ir(i,uj) be the 
unique j such that the pair (i,j) is matched under u. For each uj G fi, since 
71"^ is an internal bijection on /, it is obvious that 7r w is measure-preserving 
from the Loeb space (1,2, X) to itself. Thus, (i) of part (1) is shown. 

It is obvious that for any agent i £ I, 

(8) P ({uen:7r i (u>)=j}) = j^-j 

for any j ^ i; that is, the ith agent is matched with equal chance to other 
agents. 

Fix any i El. For any internal set C G 2q, (8) implies that Po(to G SI : vr^w) G 
C) is \C\/(N -l)i£i(£C, and (|C| - l)/(N - 1) if i G C, where |C| is the 
internal cardinality of C. This means 17 that 

(9) Po(vrr 1 (C))^^ = A (C)^A(C). 

Therefore, 7Tj is a measure-preserving mapping from (St,^, P) to (I,2q,X), 
and is measure-preserving from the Loeb space (Sl,J-,P) to the Loeb space 
(1,2, A). 18 Thus, (ii) of part (1) is shown. 
We can also obtain that 

(10) (X ®P )(tt~ 1 (C))= j P (7rr 1 (C))dX (i)^X (C)^X(C). 

A proof similar to that of Footnote 18 shows that tt is a measure-preserving 
mapping from (I x Sl,2MT, X MP) to (1,2, X). 
Next, for i^j, consider the joint event 

(11) E = {u G SI: (TTiH^-M) = (*',/)}, 



Anderson noted in [3] that (/ x Q.. T E3 T, A H P) is an extension of the usual product 
(I x ® T, A ® P). Keisler proved in [26] (see also [34]) that the Fubini property still 
holds on ( I x Q, T Kl T , A [3 P) . Thus, the Loeb product space is a Fubini extension of the 
usual product probability space. In addition, it is shown in Theorem 6.2 of [38] that when 
both A and P are atomless, (/ x Q,J H T, A IS P) is rich enough to be endowed with a 
process h whose random variables are essentially pairwise independent and can take any 
variety of distributions (and in particular the uniform distribution on [0,1]). 

17 For two hyperreals a and /3, a ~ /? means that the difference a — fi is an infinitesimal; 
see [34]. 

18 For any Loeb measurable set B £ T and for any standard positive real number e, 
there are internal sets C and D in To such that CC5CD and Aq(D — C) < e. Thus 
^(CTlCTrr^CTrr \d), and 

Po^r^P) - Trr 1 ^)) A (P - C) <£, 

which implies that tt~ 1 (B) is Loeb measurable in T. Also, A(C) < P(7r~ 1 (P)) < A(P/), 
and thus |P(7r~ (B)) — A(P)| < A(P/ — C) < e for any standard positive real number e. 
This means that P(n^ 1 (B)) = A(B). 
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that is, the ith agent is matched to the i'th agent and the jth agent is 
matched to the j'th agent. In order to show the measure-preserving property 
of the mapping (tTj, 7Tj) in the following paragraph, we need to know the value 
of Po(E) in three different cases. The first case is = (j, i), that is, the 

ith agent is matched to the jth agent and the jth agent is matched to the ith 
agent. In this case, Pq{E) = 1/(N — 1). The second case is that of Pq{E) = 0, 
which holds when i' = i or j' = j (the ith agent is matched to herself, or the 
jth agent is matched to herself), or when i! = j but j' ^ i (the ith agent is 
matched to the jth agent, but the jth agent is not matched to the ith agent), 
or when j' = i but i' ^ j (the jth agent is matched to the ith agent, but the 
ith agent is not matched to the jth agent), or when i! = j' (both the ith agent 
and the jth agent are matched to the same agent). The third case applies if 
the indices i, j and i! ,j' are completely distinct. In this third case, after the 
pairs (i,i') 5 (j, j') are drawn, there are iV — 4 agents left, and hence there are 
(N — 5)!! ways to draw the rest of the pairs in order to complete the matching. 
This means that P (E) = (N - 5)!!/(iV - 1)!! = 1/((N - 1)(N - 3)). 

Let (/ x I,I (g)2"o,Ao <8> Ao) be the internal product of (I,Io, Ao) with 
itself, and (/ x 1,1 MI, A M A) the Loeb space of the internal product. Fix 
any i, j G / with i ^ j. Let D be the diagonal :i' G I}. The third case 

of the above paragraph implies that for any internal set G G2o 

P ({u G n : (TTiH.Tr^o;)) £G-(DU ({i, j} x /) U (/ x {i, j}))}) 

|G-(DU({*,j}x/)u(Jx{i,j}))| \G\ nm . vr . 

By using the formula for Po(E) in the first two cases, we can obtain that 
P ({u G : fciu),^)) G (D U ({i, j} x /) U (/ x {i, j}))}) 

(13) 

= ~0. 

N-l 

Equations (12) and (13) imply that 

(14) P ({uj eQ:(TTi(u}),TTj(u)) G G}) ~ (A ® A )(G). 

A proof similar to that of Footnote 18 shows that (7rj,7r,-) is a measure- 
preserving mapping from (£l,J-, P) to (J x I, I MI, A M A). 

Let a be an X-measurable type function with a distribution p on S, 
Ik = a _1 ({/c}) and pk = A(I&) for 1 < k < K . Let g = a(ir). Then, for any 
1 < k < K , g -1 ({/c}) = 7r~ 1 (/fc), which is Loeb product measurable in IMJ- 
with A Kl P-measure p^ because of the measure-preserving property of ir. 
Hence, g is I M X'-measurable. For each i G I, the measure-preserving prop- 
erty of 7Tj implies that gi has the same distribution p as a. 
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Fix any i,j £ I with i 7^ j. For any 1 < k, I < K, the measure-preserving 
property of (7^ , 7Tj ) implies that 

P(gi = k,gj = I) =P(W ett--(K i (uj),Tr j (Lj)) 6 I k X Ii}) 

(15) 

= (A m A)(/ fe x /,) = p( 9i = k) ■ P( 9j = i), 

which means that the random variables gi and gj are independent. 19 Hence, 
part (2) is shown. 

Finally, take any Ai,A 2 € X, and let f(i,uj) = 1^ 2 ("7r(i, u;)) for all (i,uj) € 
I x £1. Then, / is an i.i.d. process with a common distribution on {0, 1} 
with probabilities A (^2) on {1} and 1 — A (^2) on {0}. By the exact law of 
large numbers in Theorem 5.2 of [38] or Theorem 3.5 of [40], 20 one has for 
P-almost all to £ f2, X Al {f Al )Z l ({^}) = A(A 2 ), which implies that X{A 2 ) = 
A(^4i n7rJ 1 (^ 2 ))/A(Ai). Hence, (hi) of part (1) follows. 

4.2. Proof of Theorem 2.6. Let M be any fixed unlimited hyperfinite 
natural number in *N 00 , and let I = {1, 2, ... , M} be the space of agents. 
Let Iq be the collection of all the internal subsets of /, and Ao the internal 
counting probability measure on Tq. Let (1,1, X) be the Loeb space of the 
internal probability space (I,1q,Xo). 

Let [3 be any I-measurable type function from / to S = {1, . . . , K}. We can 
find an internal type function a from I to S such that A({i £ / : a(i) ^ (3(i)}) = 0. 
Let A k = a~ 1 (k) and \A k \ = M k for l<k<K with J2k=i M k = M. Then 
X(A k ) = Pk — M k /M for 1 < k < K. Without loss of generality, we can as- 
sume that M k S *Noo. 21 For each k in {1, . . . ,K}, pick a hyperfinite natural 
number m k such that M k — m k £ *^oo an d q k ~ m k /M k , and such that 



1 When the type function a is internal, the type process is internal as well. However, 
the computations in (12)-(14) indicate that the random variables gi, i G I , are in general 
not *-independent; see, for example, [3] for the definition of *-independence. 

20 What we need in this paper is a special case of Theorem 3.5 in [40]. Let / be a 
process from (I x Q,T M T , A P) to a complete separable metric space X. Assume that 
the random variables fi, i £ I are pairwise independent. Then, for P-almost all ui £ £1, the 
sample function has the same distribution as / in the sense that for any Borel set B in 
X, A(/J 1 (B)) = (AHP)(/- 1 (B)). Fix any A el with A (A) > 0. Let f A be the restriction 
of / to A x f2, A A and A A H P the probability measures rescaled from the restrictions of 
A and A Kl P to {D £ T : L> C A} and {C eTB f:CC A x Q}, respectively. Then, for the 
case that the random variables fi,i € I have a common distribution /ionl, the sample 
function {f A )u) also has distribution /i for P-almost all ui £ O. 

21 When pk — 0, we may still need to divide some number mk by Mj, so that the ratio 
is infinitely close to a real number q fe . We can take £ *Noo with M^/M ~ 0. We then 
take an internal subset of / with Mk many points as Ak and adjust the rest Ai,l 7^ k, 
on some A-null internal sets. This will produce a new internal type function a with the 
desired properties. 
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N = J2i=i(M[ — mi) is an unlimited even hyperfinite natural number. It is 
easy to see that 



For each k in {1,2, . . . ,K}, let B k be an arbitrary internal subset of A k 
with m k elements, and let T > mk (A k ) be the collection of all such internal 
subsets. For given B k G V mk {A k ) f or k = 1, 2, . . . , K, let ^,b^...,b k be a 
(full) matching on I — Ufc=i B k produced by the method described in the 
proof of Theorem 2.4; there are (N - 1)!! = 1 X 3 X • ■ • X (N - 3) x (N - 1) 
such matchings. 

Our sample space Q is the set of all ordered tuples (Pi, P2, ■ . ■ , Bk, 
7r Bl ' B2 ' -' BK ) such that B k G V mk (A k ) for each fe = 1, . . . ,K, and n B u E 2,-,B K 
is a matching on I — (JfcLi Then , ^ ha s - Ilf=i £*) many ele- 
ments in total. Let Po be the internal counting probability measure defined 
on the collection of all the internal subsets of 0. Let (fi, .F, P) be the Loeb 
space of the internal probability space (f2,Fo,Po). Note that both (I,Iq, Ao) 
and (f2,Fo,Po) are atomless. Let (J x J-, XMP) be the Loeb space of 

the internal product (/ x £l,lo ® Fo, Ao <S> Po), which is a Fubini extension 
of the usual product probability space by Footnote 16. 

Let J represent nonmatching. Define a mapping tt from / x Q to IU {</}. 
For i G A fe and w = (B t , B 2 , . . . , B k ,tt Bi ' B2 '-' Bk ), if i G P fc , then ir(i,u) = 
J (agent i is not matched); if i ^ B k , then i £ I — \J^ =l B r , agent i is to 
be matched with agent tt Bi ' B2 '-> Bk (i), and let n(i,u>) = tt Bi > B2 >~' Bk (i). It 
is obvious that vrJ 1 ({J}) = [J^ =1 B r and that the restriction of tt u to / — 
U^Li Pr is a full matching on the set. Let g@ be the matched type process 
from / x Vt to Sl){J} under the type function /3; that is, g^(i,u) = (3(-ir(i,u})) 
with P(J) = J. 

When $^LiPr(l - ?r) = 0, we know that 2V/M ~ 0. For those i G A fc 
with p k > 0, it is clear that Pq{{oj G f2 : 7Tj(u;) = J}) = m k /M k ~ g fe = 1, and 
thus A K P(7r(i,uj) ^ J) = 0, which means that A8P(/(i,w) / J) = 0, and 
gf(uj) = J for P-almost all w G f2. Thus conditions (1) and (2) in Definition 
2.5 are satisfied trivially; that is, one has a trivial random partial matching 
that is independent in types. 

For the rest of the proof, assume that J2r=iPr(l — Qr) > 0. Let g be the 
matched type process from / x Q to SU {J}, defined by g(i,u>) = a(n(i,u)), 
where a(J) = J. Since both a and tt are internal, the fact that Iq <£> Fo is 
the internal power set on / x implies that g is Iq <8> Fo-measurable, and 
thus I M F-measurable. 

Fix an agent i G A k for some 1 < k < K. For any 1 < I < K, and for 
any B r G V mr (A r ), r = 1,2,..., K, let N Bi,b 2 ,...,b k bg thg number of full 

matchings on \J^ =l (A r — B r ) such that agent i is matched to some agent 



(16) 
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in A[ — B\. It is obvious that jV Sl ' S2 '"''' BK depends only on the numbers of 
points in the sets A r — B r , r = 1, . . . , K, which are M r — m r , r = 1, . . . , K, 
respectively. Hence, N^ 1 ' 2, '"' Bk is independent of the particular choices of 
B\, .E?2, • • • , Bk, and so can simply be denoted by Nu. Then, (9) implies that 

( 7) Nu _ Mi - mi 

[ ' (AT-l)!! ~ N ' 

It can be checked that the internal cardinality of the event {gi = 1} is 
\{u G U-.i G A k - fl fc ,7rf 1 ' Ba -" B *(u;) G (A t - B t )}\ 



(18) 



m k J \ V m r 

\rj^k 



for u = (5i,B 2> • • .,B Ki tt Bi ' B2 '-' Bk ). Hence (17) and (18) imply that 

n , ,\ M k -m k Nu H Mi -mi 

= l) = M k (jV-i)M * 0- ~ *)^V— 

(19) 

Al-rrH/MfiMj/M (1 - q k ) P i(l - q t ) 

~ [ qk) (N/M) ~ E£ift-(i-gr) • 

It is also easy to see that -Po({ w £ ^ : Si( w ) = <^}) = m k/M k ~ This means 
that for i G A k , 

P( 9i = l)= t~ K qk)Ml ~ qi) 



in types. Fix agents i,j G J with i j. For any 1 < Z , < 
£ r G 7> mT . (A), r = 1, 2, . . . , K, let N^ B2 '-' Bk be the nu 



Er=liMl - ?r) 

and that P(gi = J) = q k . Hence, the distribution condition on gi is satisfied 
for each i G I. 

We need to show that the random partial matching ir is independent 

, t < -fT, and for any 
y a - t uc unc number of full match- 

ings on \J^ =1 (A r — B r ) such that agents i and j are matched to some 
agents respectively in Ai — Bi and A t — B t . As in the case of N Bi,B2, '"' Bk , 
N it jt K is independent of the particular choices of B±,B2, . ■ . ,Bk and 
can simply be denoted by Nuj t . By taking G = (Ai — B{) x (A t — B t ), (14) 
implies that 

1 ' (JV-1)!!~ N N ' 

We first consider the case that both i and j belong to A k for some k in 
{1, . . . ,K}. It is easy to see that Po(<?j = J,gj = J) = m k {m k — l)/(M k (M k — 
1)), and hence that 

(21) P( 9i = J, 9j = J) = P( 9l = J)P( 9j = J)=ql 
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As above, it can be checked that the internal cardinality of the event {gi = I, 
9j = J} is 

\{u;en:ieA k -B k ,je B k ,ir? 1 -- B *(w) e (At - 

(22) 

for = (B 1 , . . .,B k ,tt Bi '-' Bk ). Hence (17) and (22) imply that 

m k (M k -m k ) Nu 



Po(9i = l,9j = J) 



M k (M k -l) (N-l)\\ 
Mi -mi 



(23) ^q k (l-q k ) 

_ gfc(l-gfc)pz(l-g;) 

E£=lPr(l-9r) 

which implies that 

(24) Pfo = i,^ = J) = P( 9i = l)P( 9j = J). 

Similarly, the events ( 9 i = J) and (gj = I) are independent. 

The event {gi = I, gj = t} is actually the set of all the ui = (B\, . . . ,Bk, 
7I .b 1 ,...,b k ^ guch th a t both i and j are in A k — B k , and agents i and j are 
matched to some agents in A\ — B\ and At — Bt, respectively. Thus, the 
internal cardinality of {gi = I, gj = i] is 

r^k 



< 25 > ( ™r n 



Hence (20) and (25) imply that 

(M k - m k )(M k -m k -l) N Ujt 



Po(9i = l,9j = t) 



(26) ~ (1 - q k 



Mfc(M fc -l) (JV-1)!! 
2 Mi - mi M t - m t 



N N 
(1 ~ gfc) 2 P;(l ~ Ql)Pt(l ~ qt) 
(Ef=iPr(l-^)) 2 

which implies that 

(27) P( 9i = l, 9j = t) = P( 9i = l)P( 9j = t). 

Hence the random variables gi and gj are independent. 
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For the case that i € A k and j € A n with 1 < k 7^ n < K , one can first 
observe that 

PC _/ _ n - ^ Mk ~ mfc ) m " j^jj 
n\9i-l,9j-J)- MkMn ( N _ 1 y A > 

v , , .x (M k -m k ){M n -m n ) N aj t 
n)\9i = hdj =t) = 



M fc M n (AT-1)!!- 

In this case, one can use computations similar to those of the above two 
paragraphs to show that the random variables </, and gj are independent. 
The details are omitted here. 

We have proved the result for the type function a. We still need to prove 
it for [3 [and for g@ = /3(7r)] . Fix any agent i € A k , for some 1 < k < K . 
For any internal set A € Xq, and for any B r € V mr (A r ), r = 1, 2, . . . , K, let 
N^' B2 '"' ,Bk be the number of full matchings on \J^ = x(A r — B r ) such that 
agent i is matched to some agent in A — U r =i ^r- Then, (9) implies that 

(29) N ^^Bk ^ \ A -[)? =1 B r \ 



(iV-1)!! N 
The internal event it J 1 {A) is, for a; = (B 1 ,B 2 , ■ ■ . ,B k ,tt Bi ' B2 '-> Bk ), 

(30) ^uen:i£A k -B k ,n Bl > B2 '---> BK (u)e U- \J B r )\. 

Hence (29) and (30) imply that 

, n b u b 2 ,...,b k 

P 0(^(A))= £ U K [M r \ \ A N _ ir 

Bfc£'Pm fc (A fc \{i}),B i eP m ;(A i ) for l^k Hr=l \m r ) v 

1 \A-\J? =1 B r 



E 



(31) 



< 



S fc €7' mfc ( J 4 fc \{t}) 1 B l e7' mi (A I ) for Z^fc Hr=l U r 



nf =1 CD * 



Let 



^ nf=iC) Mfc Af (N/M) 

(l-q k )Xo(A) 

(1 " ?fc) 



max 



1<*<X E^iPr(l-gr)* 

Then, for each i 6 I and any yl E Zq, P(tt~ 1 (A)) < c ■ \(A). Thus, Keisler's 
Fubini property as in [26] and [34] also implies that (A M P)(tt~ 1 (A)) < c ■ 
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\{A). Let B = {i el:a(i) + We know that \{B) = 0, {XMP)(tt~ 1 {B)) = 

P(ir^ 1 (B)) = for each i € Since g and g@ agree on / x — ir~ 1 (B), 
must be T M ^"-measurable. For each i G I, g% and g? agree on Q — tt~ 1 (B), 
and hence the relevant distribution and independence conditions are also 
satisfied by g@ . 

Remark 4.1. The sample space £1 in the proof of Theorem 2.6 depends 
on the choice of the internal type function a. In the proof of Theorem 3.1 
in Section 6 below, it will be more convenient to construct a sample space 
Q that depends only on the agent space I, and not on the type function a. 

Let 0, be the set of all internal bijections a from I to I such that for 
each i £ I, a(i) = i or a(a(i)) = i, and let be its internal power set. Let 
qk £ [0, 1] for each 1 < k < K. We can define an internal one-to-one mapping 
(p a from n to n by letting ip a (B 1 ,B 2 , . . . ,B k ,tt Bi ' B2 '-' Bk ) be the internal 
bijection a on I such that a(i) = i for i € UfcLi -Bfc and <r(i) = -r b ' l,b%, ""' Bk (i) 
for iEl — UfcLi-^fc- The mapping also induces an internal probability 
measure Pq* on (^j^b)- 

Define a mapping ff : (J x fj) — ► (/ U J) by letting Tt{i,u) = J if = 
i, and 7f(z,u)) = u)(i) if d)(i) 7^ i. Extend a so that a(J) = J and define 
g a : (/ x n) -> (5 U {J}) by letting 

g a {i,Cu) =Q(vf(i,a)),a)). 

Then it is obvious that 7f is still an independent-in-types random partial 
matching ir from (/ x Q,JMT,XM P a ) to / with q = (51, ... , g^) as the 
no-match probabilities. 

5. Generalized Fubini and Ionescu Tulcea theorems for Loeb transition 
probabilities. For the proof of Theorem 3.1 to follow in Section 6, we need 
to work with a Loeb product transition probability system for a sequence 
of internal transition probabilities, based on the Loeb space construction 
in [33]. We first prove a generalized Fubini theorem for a Loeb transition 
probability in Section 5.1. Then, a generalized Ionescu-Tulcea theorem for 
an infinite sequence of Loeb transition probabilities is shown in Section 5.2. 
Finally, a Fubini extension based on Loeb product transition probability 
system is constructed in Section 5.3. 

5.1. A generalized Fubini theorem for a Loeb transition probability. Let 
(/,2o,Ao) be a hyperfinite internal probability space with Iq the internal 
power set on a hyperfinite set /, and £1 a hyperfinite internal set with Tq its 
internal power set. Let Pq be an internal function from / to the space of hy- 
perfinite internal probability measures on (0, JF ), which is called an internal 
transition probability. For i & I, denote the hyperfinite internal probability 
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measure Po(i) by Poi- Define a hyperfinite internal probability measure tq on 
(/ x O,l ®^o) by letting r ({(i, w)}) = \ ({i})P 0i ({u}) for E/x fi. 
Let (1,1, A), (£l,Ti,Pi) and (IxfijIS^r) be the Loeb spaces correspond- 
ing respectively to (J,2o,Ao), (0,.Po)P(h) an d (-/" x ^,^o ®^b) r o)- The collec- 
tion {Pi : i E /} of Loeb measures will be called a Loeb transition probability, 
and denoted by P. The measure r will be called the Loeb product of the 
measure A and the Loeb transition probability P. We shall also denote ro 
by A 8> Po and rby XMP. 

When P 0i = Pq" for some hyperfinite internal probability measure Pq on 
(Q,To), To is simply the internal product of Ao and Pq, and r the Loeb 
product measure A KIP', where P' is the Loeb measure of Pq. A Fubini-type 
theorem for this special case was shown by Keisler in [26], which is often 
referred to as Keisler's Fubini theorem. The following theorem presents a 
Fubini-type theorem for the general case. 

Theorem 5.1. Let f be a real-valued integrable function on (I x£l,a(Io® 
!Fq),t). Then, (1) f. L is a (To) -measurable for each i E / and integrable on 
(Q, u(J-q), Pi) for \-almost all i E / ; (2) Jn fi(w) dPi(u) is integrable on 
(I,a(l ),X); (3) /, f n f t (co) dP t (u) d\(i) = f Ixn f(i,u) dr(i, u). 

Proof. Let Ti be the class of functions g from / x Q, to M.+ U {+00} 
that satisfy (1) for every % E I, gi(-) is cr(.Po)-nieasurable; (2) the integral 
dPi(uj) as a function from / to M + U {+00} is cr(2"o)-measurable; 
(3) fj Jq gt(u) dPi(ui) dX(i) = f Ixfl g(i,uj) dr(i,oj). It is obvious that Ti is 
closed under nonnegative linear combinations and monotone convergence. 

Now, we consider E E To ® J-q and g = 1e- Then, for each i E I, gi is the 
indicator function of the internal set E{ = {uj E O : (i, uS) E E}, which is J-q- 
measurable [and hence a( J r o)-measurable]. The integral f n gi(io) dPi(u>) = 
Pi(Ei) is the standard part °(Poi(Ei)). Since Poi(Ei) is Zo-measurable as a 
function on /, Pi(Ei) is thus o"(2o)-measurable as a function on I. Thus, 
the usual result on S'-integrability (see, e.g., [34], Theorem 5.3.5, page 155) 
implies that 

(32) / [ g t (uj)dPi(uj)d\(i)= [°P 0i (Ei)d\(i) 
JiJn J 1 

(33) =°J i P oi (E i )dXo(i) 

(34) =° To (E)= [ gdr. 
Thus, g E Ti, and hence Ti contains the algebra 2q ® .Pq. 
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Therefore TL is a monotone class. Then Theorem 3 on page 16 of [7] implies 
that Tt must contain all the nonnegative aiXo® To)-m.easx\iab\e functions. 22 

Since / is integrable on (I x Q,a(Xo ®To), r), so are both / + and /~. 
Now the fact that Jj J n ff (uS) dPi(uj) d\(i) < oo implies that for A-almost 
all i G J, Jq I^{uj) dPi(uj) < oo, and thus the <7(.Po)-iiieasurable function 
is integrable. Similarly, the <r(2o)-measurable function Jj-j (w) dPi(uS) is 
integrable since Jj J Q f^(u>) dPi(uS) dX(i) < oo. We have similar results for 
f~. The rest is clear. □ 

For any B G a (To), apply Theorem 5.1 to / = 1/ X _B to obtain that Pi(B) 
is <r(Xo)-measurable for each i G /. This means that P = {Pi : i G 1} is indeed 
a transition probability in the usual sense (see [25] ) . One can define its usual 
product A ® P with A by letting A ® P(P) = J I P i (E i ) dX(i) for each £ in 
the usual product cr-algebra cr(Zo) S^cr^o)- It is clear that (I x f2,cr(2o ® 
To), X M P) is an extension of (I x Q, a(2 ) (g> cr(^" ), A ® P). 

The following result extends Theorem 5.1 to integrable functions on (I x 
Jl,ZKlJF, r), which is the completion of (/ x Q,ct(Iq (g) To),t). 

Proposition 5.2. Let f be a real-valued integrable function on (I x Q, 
IMT,t). Then, (1) for X-almost all i El, fa is integrable on the Loeb space 
(Q,Ti,Pi); (2) fi(uj) dPi(u>) is integrable on (1,1, A); and (3) we have 



Proof. First, let E G X M T with r(P) = 0. Then there is a set A G 
cr(Xo (8> To) such that £Ci and r(A) = 0. By Theorem 5.1, for A-almost all 
i G /, Pi(Ai) = 0, which implies that P;(Pj) = 0. 

There is a real- valued o~(Io <S> .Po)- measura ble function g on / x such 
that the set E = {(i,ui) G I X Jl : f(i,uj) ^ g(i,co)} has r-measure zero. The 
above result implies that for A-almost all i€ I, is the same as gi. The rest 
is clear. □ 

Let (I',2 , X' ) be a hyperfinite internal probability space with I' Q the in- 
ternal power set on a hyperfinite set I'. One can define a new transition 
probability, {Aq ® Poi : i G /}. Define a hyperfinite internal probability mea- 
sure TQ on (/' X / x U,I' ®Tq ®Tq) by letting 



T 1 ({(z^^, W )}) = Ao(W)(Ao®^)({(i^^)}) = Ao(W)A ({^'}) J Po i (M) 



for (i',i,ui) G J' x J x J7. Then it is clear that Tq 1 is exactly the same as 
A'o <S> To- Let A' Kl r be the Loeb measure of A' ® To- By applying Theorem 
5.1, we can obtain the following corollary. 



22 There is a typo in Theorem 3 on page 16 of [7]; D should be an algebra (not a 
(T-algebra as stated). 
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Corollary 5.3. Let h be a real-valued integrable function on (I' x 7 x 
f2, c(Xq <g> Zo <S> .Fo) > A'Br). Then, the following results hold. 

(1) hi is a {Iq ® J- 0)- measurable for each i £ I and integrable on ($1, a(I ® 
.Po), A' El P) / or X-almost all i S 7 , where X' El Pi is the Loeb measure of 

(2) fj,^hi(i',uj)d\'MPi{i\Lj) is integrable on (7, <r(Xo), A). 

(3) // /// xn w) d(A' El Pi)(i', w) dA(i) = / J/xIxn i, w) d(A' B r)(i', 
i,u/). 

One can also view {Poi : i € 7} as a transition probability from 7' x 7 to 

0. Then the following corollary is obvious. 

Corollary 5.4. Let h be a real-valued integrable function on (7' x 7 x 
f2, <j(Zq <8> To ® X"o) , A' El r) . Then, the following results hold. 

(1) /i^j) is a {J- q)- measurable for each (i',i) € 7' x 7 and integrable on 
(^l, a (To), Pi) for A' IE] X-almost all (i',i) G 7' x 7, w/iere A' IE] A is the Loeb 
measure of X' (8> Ao . 

(2) J hui i)(u>) dPi(uj) is integrable on (I' x 7, u(2q <8>Zo), A' El A). 

(3) W /n ^,oH d ( A ' ® A)(i',t) = J}, xJxtl fc(t', i, w) d(A' El r)(i', 

1, u;). 

By applying Proposition 5.2, one can also extend the results in Corollaries 
5.3 and 5.4 to integrable functions on the Loeb space of (I' x 7 x fi,X ® 
Xo®.Fo,To). 

5.2. j4 generalized Ionescu-Tulcea theorem for a Loeb product transition 
probability system. In Section 5.1, the subscript is used to distinguish 
internal measures and algebras with their corresponding Loeb measures and 
Loeb algebras. In this subsection, we need to work with infinitely many 
internal measure spaces and the corresponding Loeb spaces. To avoid con- 
fusion with the notation, we shall use Q with subscripts or superscripts to 
represent internal measures; when their corresponding Loeb measures are 
considered, we use P to replace Q. Internal algebras are denoted by T with 
subscripts or superscripts, and their external versions by A with subscripts 
or superscripts. 

For each m > 1, let £l m be a hyperfinite set with its internal power set T m . 
We shall use fi n , O 00 and to denote Um=l n m, Um=i and Um=n ^m, 
respectively; also {o; m }^ l=1 , {u m }m=l and {io m }^ =n will be denoted respec- 
tively by uj n , oj°° and ct>^° when there is no confusion. 

For each n> 1, let Q n be an internal transition probability from fi™" 1 to 
(n n ,Fn); that is, for each uj n ~ l £ ft"" 1 , Q n {uj n ~ l ) (also denoted by Qf~') 
is a hyperfinite internal probability measure on (Q n ,J- n ). In particular, Q\ 
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is simply a hyperfinite internal probability measure on (£li,J-i), and Qi an 
internal transition probability from fii to (f^j-^)- Thus, Q\ <g> Q2 defines 
an internal probability measure on (Qi x 0,2, T\ <S)J~2)- By induction, Q\ ® 
Q2 <8> • • • <8> Qn defines an internal probability measure on (f2 n , <S>m=i ^m)- 
Denote Qi <8> Q 2 ® ■ • ■ <8> Qn by Q n , and (g)™ =1 JT m by JP\ Then Q n is the 
internal product of the internal transition probability Q n with the internal 
probability measure Q n ~ l . Let P n and P^u"- 1 ) (also denoted by ) be 

the corresponding Loeb measures, which are defined respectively on a{J- n ) 
and a{J- n ). Using the notation in Section 5.1, P n is the Loeb product P± £3 
P2M ■ ■ - M P n of the Loeb transition probabilities P\ , P2 , ■ . • , P n ■ 
Theorem 5.1 implies that for any set E G a(J- n ), 

(35) P n (E) = [ Pf-'iE^dP"- 1 ^- 1 ). 

That is, P n is the product of the transition probability P n with the prob- 
ability measure P n ~ l . Thus, we can also denote P n by P n ~ l M P n , and 
furthermore by K1J^ =1 P m . 

Let 7r" _1 be the projection mapping from Q n to f2 n_1 ; that is, 7r™ _1 (a;i , . . . , 
w n ) = (u>i, . . . ,uj n -i). Let F be any subset of f2 n_1 and E = F x Q n . Then, 
E^n-i = n n when uj n ~ l G F, and E^n-i = when a;™" 1 <£F.E Ee a(T n ), 
then Theorem 5.1 implies that P^ n (E^n-i) = If is ^(.F"' -1 ) -measurable 
[i.e., F G a{F n ' 1 )], and P n (£) = P n ~ l (F). On the other hand, if F G 
a(J zrn ~ 1 ), then it is obvious that E G a{T n ) and P n (£) = P n ~ 1 (F). This 
means that the measure space (Q, n ~ l ,a(J- n ~ l ),P n ~ l ) is the projection of 
(Q n , a{J rn ), P n ) under vr™ _1 . Similarly, let tr^ be the projection mapping 
from Q n to Q k for some k < n; then (O fc , a{J- k ), P k ) is the projection of 
(n n ,(j(r),P n ) under vr£. 

For a collection T> of sets and a set F, we use P x F to denote the collection 
{I? x F : D G £>} of sets when there is no confusion. Thus, cr(T n ) x 
denotes {E n x : £ n G a^)}. Let 5 = \J™ =1 [a(T n ) x which is 

an algebra of sets in One can define a measure P°° on this algebra 
by letting P°°(E n x ft£° +1 ) = P n (E n ) for each E n ea{T n ). The projection 
property stated in the above paragraph implies that P°° is well defined. Let 

The point is how to extend P°° to a countably additive probability mea- 
sure on the ex-algebra a(J-°°). This is possible by using a proof similar to that 



In fact, for each uj n £ O™, we have 

n n 

Qm({(wi,. ■ .,Wn)}) = Qm(Wl,. . . , W m _i) ({«„}). 
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of Proposition 3.3 of [39]. The result is a version of the Ionescu-Tulcea the- 
orem (see [25], page 93) for the Loeb product transition probability system 

{p 1 mp 2 m---®p n }™ =1 . 

Theorem 5.5. There is a unique countably additive probability measure 
on a(J-°°) that extends the set function P°° on £; such a unique extension 
is still denoted by P°° and by M^ =1 P m . 

Proof. Let {C n }^ =1 be a decreasing sequence of sets in T 00 with empty 
intersection. By the construction of J 700 , one can find a sequence of internal 
sets {A n }^ =1 and a nondecreasing sequence {k n }^ =1 of nonnegative integers 
such that C n = A n x and A n G T kn . For £ < n, let n^. 1 be the mapping 

from Q kn to Q kl by projecting a tuple in Q fcn to its first ki coordinates; 
then ir* l n (A n )CA e because {C n }^ = i is a decreasing sequence of sets. Take 
the transfer {A: n } ne *pj of the sequence {k n }™ =1 , and the respective internal 
extensions {A n } n <z*^ and {il n } ne »N of the internal sequences {A n }^ =1 and 
{r2 n } ng N. By spillover and ^i-saturation (see [34]), one can obtain h G *Noo 
such that for all n<h, A n C O fe ™ and vr^ (A n ) C A e for all £ G *N with I < n, 

where 7T^ ! is defined in exactly the same way as in the case of finite n. 

We claim that A n = for all n G *Noo with n < h; if not, one can find 
such an n with u) kn = (ui,... ,u) n ) G A n . Then uo kl G A\ for any i G N. If 
k n G *Noo, then it is obvious that {uj m }^ =1 is in Q for all I G N, which 
contradicts the assumption that the intersection of all the Ce is empty. If 
k n G N, one can choose io m arbitrarily for any m > k n to obtain the same 
contradiction. Hence the claim is proven. 

By spillover, we know that for some n G N, A n = 0, and so is C n . Thus, we 
obtain a trivial limit, limj^oo P°°(C n ) =0. This means that P°° is indeed 
countably additive on T°° '. As in [33], the Caratheodory extension theorem 
implies that P°° can be extended to the a-algebra o~{J-°°) generated by J-°° , 
and we are done. □ 

The following result considers sectional measur ability for sets in o(J-°°). 

Proposition 5.6. Let G be a a (JF 00 ) -measurable subset of Q°° . Then, 
for any {w m }^ =1 G , the set G^ +i = {J n G Q n : (uj[, . . . ,u' n ,u n+1 , 
uJn+2, . . .) G G} belongs to a(T n ), while the set 

G w n = G fij^ '■ • • • > w n,^ n+ i,^ n +2' ■ ■ ■) £ G} 

belongs to o-(U™ =n+1 [(®T =n+1 F k ) x ^C+iD- 

Proof. The collection of those sets G in a{J-°°) with the properties is 
clearly a monotone class of sets and contains the algebra J-°°; and hence it 
is a {J- 00 ) itself by Theorem 1 on page 7 of [7]. □ 
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The following corollary follows from Proposition 5.6 immediately. 

Corollary 5.7. Let ix n be the projection mapping from fi 00 to [i.e., 
7^(0;°°) =6t) nj_ Then the measure space (fi n , a(J- n ), P n ) is the projection of 
the measure space (f2°°, o~(J-°°), P°°) under ir n in the sense that for any 
F C n n , F € a(T n ) if and only if (tt 71 )" 1 ^) = F x G c^F 00 ) mi/i 

P n (F) =P 00 ((vr n )- 1 (F)). 

5.3. Fubini extensions based on a Loeb product transition probability sys- 
tem. Let (I,Xo,Ao) be a hyperfinite internal probability space with Iq the 
internal power set on a hyperfinite set I. We follow the notation and con- 
struction in Section 5.2. Denote (f2o> F"o, Qo) = (F^O)Ao) an d repeat the 
process of constructing a countably additive measure K^ =0 P m on (J x 

The following lemma is a restatement of Keisler's Fubini theorem to the 
particular setting. Since the marginals of K]^ =0 P m on L and Vl n are respec- 
tively A and P n , we can write Kl^—pPm as A KIP™. 

Lemma 5.8. For any n>l, the space (I x J7 n , a(T ®J rn ) 1 Kl™ =0 P m ) is a 
Fubini extension over the usual product of (I, c(Tq), A) and (Q. n ,a(J- n ),P n ) . 

The following is a Fubini- type result for the infinite product. 

Proposition 5.9. The space (I x ,a(\J^ =l (T Q ® T n ) x n%> +1 ), 
M^ =0 P m ) is a Fubini extension over the usual product of the probability 
spaces (1,1, A) and (ft 00 ,^ 00 ),^ 00 ). 

Proof. We only check that the Fubini-type property holds for any set 
E G {J^ = i(Tq ® J 7 ™) x fiJJ+i- The rest of the proof is essentially the same as 
that of Theorem 5.1. 

It is clear that there exists a set FgIo® J- n such that E = F x f^J^i- 
By definition, K~ =0 P m (F) = =0 P m (F). By Lemma 5.8, 

1 P m (F) = \®P n (F)= f P n (F i )dX(i)= f \(F u n)dP n (u n ). 

m=0 J J jQn 

On the other hand, Fj = Fj x and F^oo = F^n . By the fact that 

F°°(Fj) = P n (Fi), the projection property in Corollary 5.7 implies that 

i F m (F) = f P°°(F i )dX(i) = f X(F w o )dP 00 (u; 00 ). 

"1=0 J J Jqoo 

This means that the Fubini property does hold for sets in U^i(^o <S»F n ) x 
00 n 
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For simplicity, we denote (j(Iq), a(T n ), o(Iq® J-'n), fJ (^ rn ) 5 ^(Iq ® J^" 1 ), 
a(F°°) and a({J™ =l {l ® ^ n ) x ft~ +1 ), respectively by I, A n , IE An, A n , 
IMA n , A°° andl^A 00 . 

We restate some of the above results using the new notation. Corollary 
5.7 implies that (n n ,A n ,P n ) and (I x tt n ,IMA n , XMP n ) are the respective 
projections of (£1°° , A°° , P°°) and (I x ft 00 , 1 El „4°°, A El P°°). Since (I x S7 n , 
IMA n , A Kl P n ) is the Loeb product of two Loeb probability spaces (1,1, A) 
and (£l n , A n , P n ) , it is a Fubini extension of the usual product probability 
space. In addition, the Fubini property in Proposition 5.9 says that (I x 
Q°°,I M A°°, A M P°°) is a Fubini extension of the usual product of the 
probability spaces (1,1, X) and (ft 00 , A°°, P°°). 

6. Proof of Theorem 3.1. Let (p°, b, q, v) be the given parameters for the 
dynamical system B. Let M be a fixed unlimited hyperfmite natural number 
in *Noo, I = {1, 2, . . . , M}, Iq the internal power set on /, and Ao the internal 
counting probability measure on Iq. Let a : 1 — ► S = {1, 2, . . . , K} be an in- 
ternal initial type function such that Ao(a° = k) ~ pk for each k = 1,. . . , K. 
What we need to do is to construct a sequence of internal transition prob- 
abilities and a sequence of internal type functions. The results in Sections 
5.2 and 5.3 can then be applied to obtain a Loeb product transition proba- 
bility system. Since we need to consider random mutation, random partial 
matching and random type changing at each time period, three internal 
measurable spaces with internal transition probabilities will be constructed 
at each time period. 

Adopt the notation used in Section 5.2. Suppose that the construction 
for the dynamical system B has been done up to time period n — 1. Thus, 
{(fl m , Tm, Qm)}m=i and {a'jfJo 1 nave been constructed, where each ft m is 
a hyperfmite internal set with its internal power set J~m, Qm an internal 
transition probability from ft™" 1 to (£l m ,F m ), and a 1 an internal function 
from / x Q 51 to the type space S. 

We shall now consider the constructions for time re. We first work with 
the random mutation step. For each 1 < k < K, p^ is a distribution on S 
with pk(l) =bki, the probability for a type-fc agent to mutate to a type-/ 
agent. Let = S (the space of all internal functions from I to S) with 

its internal power set Tzn-2- 

For each i £ I, u 3n ~ 3 G ft 3 " -3 , let 7^ 3n ~ 3 = p a -n.-iu^n-zy That is, if 
a n_1 (i,u; 3n ~ 3 ) = k, then jf 3 " 3 = p^- Define an internal probability measure 
Q3n.-2 011 (S 1 ,Fzn-2) to be the internal product measure ILe/ 7f " - 25 Let 
h n ■ (I x ]lm=i 2 ^m) ^Sbe such that /i n (i,o; 3n ~ 2 ) = LU 3n - 2 (i). 



This is possible since A is atomless. 
25 A hyperfmite product space is a common construction in nonstandard analysis, whose 
coordinate functions also give a hyperfmite sequence of *-independent random variables. 
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Next, we consider the step of random partial matching. Let (f23 n _x, J~3n-i) 
be the internal sample measurable space in Remark 4.1. 

For any given u; 3n_2 G £l 3n ~ 2 , the type function is h^ jin _ 2 {-) , denoted by 

3n — 2 

a for short. Let Q% n _i be the internal probability measure correspond- 
ing to the internal probability measure Pff in Remark 4.1. Define a map- 
ping ir n : (I x ft 3 ™- 1 ) — > (I U J) by letting vr n (i, cj 3 ™" 1 ) = vf (i, w 3 n-i); thus, 
% n (i,ij Zn ' x ) = J if w 3n _i(i) =i, and 7r n (z,a; 3ri-1 ) =w 3n _i(i) if w 3n _i(z) 7^ i. 
Extend h n so that /i n (J, u; 3n ~ 2 ) = J for any oj 3n ~ 2 G fi 3n " 2 . Define 5™ : (J x 
SI 3 ™" 1 ) -> (S U {J}) by letting 

g n (t, u; 3 ™" 1 ) = /i n (vr n (i, u; 3 ™- 1 ) , u 3n - 2 ) , 

which means that g n (i, w 3n_1 ) = g Q (z,u;3 n „i). 

Finally, we consider the step of random type changing for matched agents. 
Let Q,3 n = S with its internal power set J~s n ; each point u>s n € is an 
internal function from / to 5. For any given u; 3n_1 G il 3n_1 , the space / of 
agents is divided into K 2 + K classes: those in type k who are not matched, 
or matched to some type-/ agents. For 1 < k, I < K, v\~\ is a distribution on 
S and Vklir) the probability for a type- A; agent to change to a type-r agent 
when the type- A: agent meets a type-/ agent. 

Define a new type function a n :(I x Sl 3n ) — ► S by letting a n (i,uj 3n ) = 

Fix uj 3n ~ l G il 3n_1 . For each i El, (1) if a> 3n _i(z) = i (i is not matched at 
time n), let be the probability measure on the type space S that gives 

probability 1 to the type h n (i, ui 3n ~ 2 ) and zero for the rest; (2) if w 3n _i(i) ^ i 
(i is matched at time n), h n (i,uj 3n ~ 2 ) = k and h n (ui3 n -i(i),uj 3n ~ 2 ) = /, let 
be the distribution v^i on S. Define an internal probability measure 

Q%n 1 011 ^ to be the internal product measure Ylizi T t " 1 • 

By induction, we can construct a sequence {(f2 m , ,F m , Q m )}™ =1 of internal 
transition probabilities and a sequence {a °f type functions. By using 
the constructions in Sections 5.2 and 5.3 via an infinite product of Loeb 
transition probabilities, we can obtain a corresponding probability space 
(J x O^IB^.ABP 00 ). 

From now on, we shall also use (£l,J-,P) and (/ x J-, A M P) to 

denote (tl 00 , A°° , P°°) and (I xQ 00 ,imA°° , XMP°°), respectively. Note that 
all the functions, h n ,TT n ,g n ,a n , for n = 1,2, . . . , can be viewed as functions 
on I x Q, and h n ,g n ,a n are 2 M .^-measurable for each n > 1. 



In general, for a hyperfinite collection {Xi}f =1 (with n infinitely large) of *-independent 
random variables on an internal probability space (Q,^, Po) with mean zero and variances 
bounded by a common standard positive number C, the elementary Chebyshev's inequality 
says that for any positive hyperreal number e, Po(|Ai + ■ ■ ■ + X n \/n > e) < C/ne 2 . By 
taking e = l/n 1//4 , it implies that for nearly all weft, |Ai + • • • + A„|/n ~ 0, which was 
also noted in [26], page 56. 
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We still need to check that our internal constructions above lead to a dy- 
namical system D with random mutation, partial matching and type chang- 
ing that is Markov conditionally independent in types. We assume that the 
conditions for random mutation, partial matching and type changing as well 
as Markov conditional independence in types are satisfied up to time n — 1. 
As in the proof of Lemma 6 in [15], Lemma 5 in [15] implies that the random 
variables a™ -1 and a™ -1 are independent for i ^ j. It remains to check the 
conditions for random mutation, partial matching and type changing as well 
as Markov conditional independence in types for time n. 

For the step of random mutation at time period n, we have for each agent 
i G /, and k,l G S, 

p(/ i r=/,«r l =A ; ) 

= P 3n ~ 2 (A 3n . 2 ) 

(36) = p3n-2(| u; 3n-3 £ Q 3n-3 . = fc | x ^7-{i} x {/}«)) 

p k (l) dP 3n - 3 (oo 3n - 3 ) = buPiar 1 = k), 



'{«™-V 3 ™- 3 )=fc} 

where 

h n (i,UJ 3n -' 3 ,UJ3n-2) =W 3n _ 2 (i) =/}, 

which implies that (1) is satisfied. 

When j El, it is obvious that for any l,r £ S, and any (a°, . . . , a™ -1 
and (a*-, . . . , a™ -1 ) in S n , for the event 



B 1 = {h2 = l,K = r,(al...,ar 1 ) = (al...,a 



n—l\ 

j ' i v^i i ■ ■ ■ i ; v""i ' • ■ • j "i /j 

a°,...,^- 1 ) = (a J ,...,a J "- 1 )}, 



we have 



P{B l ) = P 3n ~ 2 {B 2 



>2) 

(37) = P 3n ~ 2 {B 3 x (/ x S J -W} x x {r}W)) 

= P 3 "- 3 ( J B 4 )-p ar i(/)-p arl (r), 

where 

B 2 = {( W 3 ™- 3 , W3n „ 2 ) G n 3n ~ 2 : (a?, . . . , <" ^(a; 3 "" 3 ) = (a°, . . . , a?" 1 ), 

(a°,...,a J "- 1 )( W 3 ™- 3 ) = (a J °,...,a^ 1 ), 

^3n-2(«) = l,U3n-2(j) = r}, 

B 3 = {u 3n ~ 3 G n 3n ~ 3 :(<*?,..., ar 1 )^ 3 "" 3 ) = (a°, . . . , a?" 1 ), 

(a J °,...,a"- 1 )(a; 3 "- 3 ) = (a J ,...,^- 1 )}, 
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and 

B 4 = {a, 3 "" 3 G n 3n ~ 3 : (a° m , a^ 1 )^ 3 ) = (a° m , a^ 1 ), m = 

Equations (36) and (37) imply that for any l,r £ S, 

P(h? = I, h] = r\(a° m , a™" 1 ), m = 

= P(h? = l\ar 1 )P(h] = r\a]- 1 ). 

Hence (5) in the definition of Markov conditional independence for random 
mutation is satisfied. 

Equation (37) together with the independence of a™ -1 and aij~ implies 

that /i"(-) and hj(-) are independent. As in (2), the expected cross- 

sectional type distribution immediately after random mutation. The exact 
law of large numbers in [38] and [40] (see Footnote 20 above) implies that 
for P 3n - 2 -almost all 0J 3n ~ 2 G 3n ~ 2 , 

(38) A ({,'G/:/ i ^- 2 (o=n)=pr 1/2 

for any I G 5. 

For the step of random partial matching, the definition of ir n clearly 
shows that for each ui G f2, the restriction of vr™(-) to / — (tt™)~ 1 ({J}) is 
a full matching on that set. We need to check that the function g n : (I x 
il) — > (/ U {J}) satisfies the required distribution and Markov conditional 
independence conditions. 

For any given w 3n ~ 2 G VL 3n ~ 2 , take a = h n 3n _ 2 (-) as in the construction 
Qzn-i above. Then, for each i G / and c G 5 U {J}, 

(39) 

= P^({Oen:g a (i,0)=c}). 
Moreover, for each j G I with j ^ i and d G S U {J}, 

Q3n-i 2 ({^3n-i G 3n -i : a, 3 "" 1 ) = c,g n (j,u; 3n ~ 1 ) = d}) 
(40) = P a ({cD G : g a (i, w) = c, g a (j, S>) = d}) 

~^(^ = C ).fg*(^ = d). 
For each agent i € I, k,l & S, (39) implies that 
P(tf = J,h? = k) 

= P 3 *" 1 ^ 3 "- 1 G ft 3 '"" 1 : /i n (i,cj 3n " 2 ) = jfe.^^w 3 "" 1 ) = J}) 

(41) 

= / gfcdP 3 "" 2 ^- 2 ) 

./{/i»(i,cj 3 "- 2 )=fc} 
= q k P(h? = k). 
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In addition, we obtain from (38) that 
P(gf = l,hf = k) 

= P 3n - 1 ({u 3n - 1 € ft 3 "- 1 : h n (i,uj 3n - 2 ) = ife,0 n (i > w 3n ~ 1 ) = I}) 



{h n (i,u} 3n - 2 )=k} 

(i- qk )(i- qi )X({i'el:h^ 2 ( i ') = l}) dp3n _ 2 2> 

{/l » (iiW s»- a)=fc} E f =1 (i _ ?r ) A ({z' g I : fc» (*/) = r )} 



(1 - q k ){l - qi)p" 1/2 



p(K = k) 



(i-%)(i-%)pr 1/2 



/.. s_n-l/2 
Er^lC 1 -9r)Pr 

Hence, (41) and (42) imply that (3) holds. 

Fix i^j gl. Take any (a?, . . . , a™" 1 ), (o§, . . . , a™" 1 ) g S n , l,r E S, and 
c, d 6 5 U { J}. Let D be the set of all uj 3n ~ 2 € 3n " 2 such that /i™(cj 3n - 2 ) = /, 



a™" 1 )^ 3 ™- 3 ) = (a°, . . . ,a]~ 1 ). Then, (38) and (40) imply that 
P(g? = c,<$ = d,h? = l,h] = r, 

(43) = / Pg^^OdP^-V^ 2 ) 

Jo 

= / P h Z^-i(D 2 )-P h Z^-i(D 3 )dP 3n - 2 (oj 3n - 2 ) 
Jd 

= P 3n ~ 2 {D)P{g™ = c\h? = l)P{tft = d\W = r), 



where 



D! = {wsn-i G flan-l : <?"(*, c^" 1 ) = c, 5 "(j, w 3 "" 1 ) = d}, 
Z) 2 = € : 3 « 3 »- 2 (i, tD) = c}, 

D 3 = {cD G : A 3 - 2 (j» = d}, 

which means that the Markov conditional independence condition as formu- 
lated in (6) for random partial matching is satisfied. 

Finally, we consider the step of random type changing for matched agents 
at time n. For k g S, let v k j be the Dirac measure at k on S, that is, 
u kJ (r) = 5% for each r g S. If h n (i, oj 3n ~ 2 ) = k for k g S and ^(w 3 ™- 1 ) = c 
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for c G S U { J}, then the measure rf 3 " 1 in the definition of Q 3 ^ 1 on S 1 
is simply v kc . 

For each agent i El, and for any r, A; G S, and c € S 1 U {</}, we have 

P(o? = rA ? = fe,S? = c) 

= P 3n ({(u 3n -\ u 3n ) G ft 3n : h?(uj 3n - 2 ) = k, 
(44) 9 ™(^ 1 ) = C , W3 n«=r}) 

= / u kc {r)dP 3n -\u 3n - 1 ) 

./{u;3n-l e fi3n-l . ^ I (w 3 ™- 2 )=fc,g^(^ 3 "- 1 )=c} 

= Z / fcc (r)P(^ = fc, 5 r = c), 

which implies that (4) is satisfied. 

Fix i^jel. Take any (a°, . . . , a? -1 ), (a°, . . . , a™" 1 ) G 5 n , k,l,r,t G 5, and 
c,d€: SU {J}. Let £ be the set of all cj 3 ™" 1 G ft 3 "- 1 such that ^(a; 3 ™" 1 ) = 
c, ^(w 3 " -1 ) = d, ^(o; 3n - 2 ) = k, h]{u 3n ~ 2 ) = I, (aj,...,a^ 1 )(w 3n - 3 ) = 
(a?, . . .,<-") and (a$, . . . ,a n ~ 1 ){uj 3n - 3 ) = (a$, . . . ,a n ~ 1 ). Then, letting 

El = {a? = r,a? = t,g? = c,g] = d,h? = k,h] = I, 

(a m ,...,a m ) (a m ,...,a m ),m 

(44) implies that 

P(E 1 )= [ P^\W3nett 3n :u; 3n (i) = r,Lu 3n (j)=t})dP 3n - 1 (u; 3n - 1 ) 

JE 

(45) = f v kc {r)m{t)dP 3n -\u 3n - 1 ) 

JE 

= P(E)P(a? = r\h? = k, gf = c)P(a] = t\h] = I, g] = d), 

which means that the Markov conditional independence condition as formu- 
lated in (7) for match-induced random type changing is satisfied. 

Therefore, we have shown that D is a dynamical system with random 
mutation, partial matching and type changing that is Markov conditionally 
independent in types. 
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